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ABSTRACT 


Partitioned  solution  procedures  for  direct  time  integration  of  second- 
order  coupled-field  systems  are  studied  from  the  standpoint  of  accuracy.  Such 
procedures  are  derived  by  three  formulation  steps:  implicit  integration  of 
coupled  governing  equations,  algebraic  partitioning,  and  extrapolation  to  ad¬ 
vance  the  solution.  It  is  shown  that  the  combined  effect  of  partition,  extra¬ 
polation  and  computational  paths,  dictates  the  choice  of  stable  extrapolators 
and  the  preservation  of  rigid-body  motions.  Stable  extrapolators  for  various 
conputational  paths  are  derived  and  implementation-extrapolator  combinations 
which  preserve  rigid-body  motions  are  identified.  A  spectral  analysis  shows 
that  the  partition  Itself,  when  stable,  affects  primarily  frequency  distor¬ 
tion.  Finally,  the  advantages  and  shortcomings  of  five  partitioned  procedures 
are  discussed  as  a  guide  to  practical  applications. 
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INTRODUCTION 


The  computer  simulation  of  coupled-field  dynamic  systems  generally  In¬ 
volves  two  distinct  phases.  A  spatially  discretized  model  Is  produced  through 
finite  element,  boundary  element  or  finite  difference  techniques.  The  re¬ 
sulting  equations  of  motion  are  then  solved  by  direct  time  Integration. 

Often  the  response  characteristics  of  the  component  (single-field)  sub¬ 
systems  are  markedly  different.  As  examples  of  this  situation  we  can  cite 
fluid-structure,  soil-structure,  and  certain  structure-structure  Interaction 
problems.  This  difference  may  be  due  to  material  properties,  different 
spatial  discretization  techniques,  localized  nonlinearities,  boundary  layer 
effects,  or  other  computational  considerations. 

For  coupled  systems  befitting  the  preceding  description.  It  Is  natural 
to  think  of  tailoring  the  time  integration  procedure  to  subsystem  response 
features.  An  efficient  way  of  Implementing  this  strategy  Is  to  view  each 
single-field  subsystem  as  a  computational  entity,  which  Is  treated  by  an  ap¬ 
propriate  time  Integration  program  module  ("analyzer").  The  time-advancing 
process  results  from  coordinating  the  execution  of  the  analyzers  in  sequen¬ 
tial  or  parallel  fashion.  Field  coupling  effects  are  Incorporated  through 
cyclic  Information  transfers  predicated  throuqh  appropriate  extrapolation 
formulas.  Solution  procedures  based  on  this  approach  were  called  partitioned 
transient  analysis  procedures  (or  simply  partitioned  procedures)  In  [1],  In 
this  paper,  a  general  partitioning  procedure  was  presented  for  treating  second 
order  coupled-field  equations  of  motion. 
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The  rational  selection  of  a  specific  partitioning  procedure  among  com¬ 
peting  ones  demands  an  understanding  of  numerical  stability  and  accuracy 
characteristics.'  In  [1],  a  general  stability  analysis  technique  was  developed 
to  eliminate  case-by-case  stability  considerations.  This  technique  has  pro¬ 
vided  a  coumon  basis  for  evaluating  procedures  on  grounds  of  stability,  which 
Is  critical  in  the  Initial  method-design  stages. 

Once  satisfactory  stability  requirements  are  achieved,  it  remains  to 
assess  accuracy.  The  understanding  of  this  subject  has  so  far  been  based 
largely  on  numerical  experiments  [2-7].  This  paper  endeavors  to  fill  the 
void  by  developing  an  accuracy  assessment  theory.  More  specifically,  the 
following  three  points  are  addressed:  (a)  strong  interdependency  of  accuracy, 
extrapolator  formula,  and  computational  sequence  details;  (b)  preservation  of 
rigid-body  motions  under  partitioning;  and  (c)  Incremental  effect  of  parti¬ 
tioning  on  numerical  damping  and  frequency  distortion  measures.  These  three 
aspects  provide  sufficient  grounds  for  a  comparative  assessment  of  competing 
stable  procedures. 
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REVIEW  OF  PARTITIONED  PROCEDURES 


The  three  basic  formulation  steps  discussed  In  [1]  are:  Implicit  time 
discretization  of  the  entire  system  of  coupled  field  equations;  partitioning 
of  the  coefficient  matrix  (herein  termed  solution  matrix)  that  premultiplies 
the  current  state  vector;  and  application  of  extrapolators  to  the  current- 
state  terms  that  end  up  on  the  right-hand  side  of  the  partitioned  algebraic 
equations.  The  resulting  partitioned  procedure  embodies  those  previously 
proposed  by  Belytschko  and  Mullen  [2,3],  Park,  Fellppa  and  DeRuntz  [4],  and 
Hughes  and  Liu  [5,6].  The  procedure  Is  amenable  to  a  root-amplification 
stability  analysis  that  eliminates  the  need  for  ad-hoc  Investigation  of 
special  partitions.  Below  we  Surmarlze  the  equations  from  [1]  that  will  be 
needed  for  the  ensuing  exposition. 


Equations  of  Motion: 


M  ii  ♦  0  i  ♦  K  u  ■  f 


Implicit  Integrator: 


4  -n 


m 


?n-J> 


Impl 1 c 1 t 
Difference 
Equations  : 


E  u  •  q 
~  —  n 

E  •  M  ♦  6D  +  62  K 

-Srf  A#  A»  0^ 

9  n  *  *2  f  „  4  M  (*' ♦  «  »>“  ) 
J  n  — n  —  —  n  — n 


(1) 


(2) 


(3) 
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Partitioned 

Difference 

Equations: 


£  •  ai  * 

<■  £.1  *  £2 

E  •  Ei  *  h 

£1  •  S“£i  *  £., 

£2  ■  4  £2  ♦  >2  £2 

E  1  u  *  g  -  E  -»  u^ 
~l-n  <*  n  -2  -n 


(4) 


In  the  foregoing  equations,  M,  D  and  K  are  mass,  damping  and  stiff¬ 
ness  matrices,  respectively;  u  and  f  are  solution  and  appl led- force  N-vectors, 
respectively;  dot  superscript  denote  temporal  differentiation;  a  is  a  gener¬ 
alized  time  stepsize,  a,  and  e<  are  coefficients  of  the  linear  multistep  in- 

J  J  /q\ 

tegration  formula  used,  and  hu  is  the  corresponding  historical  vector;  uVHi 

denotes  predicted  or  extrapolated  solution  vector.  Subscript  n  indicates  com¬ 
puted  values  pertaining  to  the  sample  time  t  ■  tR.  Finally,  subscripts  1  and 
2  identify  terms  associated  with  the  sum-partition  of  matrices  D  and  K. 


Table  1  provides  a  set  of  examples  of  specific  partitions  of  the  stiffness 
matrix  K_  tailored  to  a  coupled  mechanical  system  consisting  of  two  fields: 
x  and  y,  which  interact  through  a  boundary  b. 
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EFFECT  OF  COMPUTATIONAL  PATH  AND  EXTRAPOLATOR  SELECTION 


Ref.  [8]  discussed  the  organization  of  the  time-advancing  calculations 
and  identified  four  "computational  paths",  which  were  labeled  (O'),  (0),  (1) 
and  (2).  Paths  differ  according  to  the  way  in  which  auxiliary  quantities  such 
as  velocities  are  calculated.  In  [1]  it  was  shown  that  computational  paths 
can  have  a  significant  effect  on  the  stability  of  partitioned  solution  pro¬ 
cedures,  and  that  this  effect  is  tied  up  with  the  selection  of  extrapolators . 
For  example,  when  the  trapezoidal  Integration  rule  is  used  in  conjunction  with 
the  last-solution  extrapolator , 


the  solution  remains  stable  if  the  velocity  is  computed  by  integrating  the 
acceleration  that  is  in  turn  obtained  from  the  equations  of  motion.  On 
the  other  hand,  if  the  velocity  is  computed  from  the  difference 
formula  (2),  the  solution  becomes  unstable.  In  terms  of  the  nomenclature  of 
[8],  these  two  cases  correspond  to  the  ( CO  * )  and  the  (Cl ^computational  path, 
respectively,  as  sumrvsrized  in  Table  2. 


Combinations  of  different  extrapolators  and  computational  paths  will  give 
rise  in  general  to  different  characteristic  equations,  even  if  all  other 
factors  (integration  formula,  partition)  remain  the  same.  This  has  major  im¬ 
plications  as  regards  stability  and  accuracy.  To  study  these  effects,  we  in¬ 
troduce  a  fairly  general  class  of  extrapolators: 


(<* i  u 


J  -"-J 


*66 


J  ^-J 


(6) 
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where  <» j ,  and  are  nunerlc  coefficients.  This  class  Is  broader 
than  that  considered  In  [1],  which  did  not  account  for  historical  derivative 

The  characteristic  matrix  equations  for  the  four  computational  sequences 
of  Table  2  can  be  derived  by  seeking  nontrivial  solutions  of  the  form 


xk  u. 


(7) 


The  resulting  formulas  are  collected  In  Table  3.  An  undamped  system 
(0  •  0)  Is  assumed  for  simplicity;  the  damped  case  can  be  handled  In  a 
similar  fashion,  for  derivation  details  Appendix  A  of  [1]  should  be  con¬ 
sulted. 

Nunerlcal  stability  requires 

Re(Zj)  <  0  ,  1  •  1 ....  2  Mm  (8) 

for  each  characteristic  equation  In  Table  3.  This  Imposes  In  turn  con¬ 
straints  on  the  selection  of  extrapolators  for  each  path;  re contended  ex- 
trapolators  are  listed  In  Table  2.  For  the  Important  case  of  the  trapezoidal 
Integration  rule,  stable  extrapolators  are  listed  In  Table  4. 

Remark:  If  partitioning  Is  performed  at  the  differential  equation  level  as 
opposed  to  the  present  difference  equation  level,  viz., 

*  «  4  $  ♦  £l  M  *  l  *  £2  •  L?  “  (9> 

the  following  characteristic  equation  results  with  0  ■  0 

det  I  o2  H  ♦  62  ♦  /  (o?  (,  ♦  o  o  e,  ♦  p?  «2)  I  •  0  (10) 
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TABLE  4.  STABLE  EXTRAPOLATES  FOR  DIFFERENT  COMPUTATIONAL  PATHS 
FOR  TRAPEZOIDAL  INTEGRATION  FORMULA 


Computational 

Path  ; 

Stable  Extrapolator 

(CO*) 

(Vi* 

U(p)  -  ) 

"  )  h2  ~  f 

l<Vi  * h  Vi 4  rVi' 

(Cl) 

„<*> 

„  ■  Vi  *  h  Vi 

(C2) 

2 

U(p)  .  un  l  ♦  h  0n.,  ♦  h,  V, 

n 

(JO) 

u(p)  .  h 
n  Vi  4  ?  Vi 

•  This  Is  used  extensively  In  [1]. 

♦  Applicable  strictly  to  Implicit-explicit  partitions  with  somewhat 

reduced  stability  limits,  e.g.,  «  •  h  <  1.57  for  DOF-by-DOF 

I  -  E  partitions. 

Integrator:  ■  ufll  ♦  6  (0n  ♦  ).  6  •  j 

Extrapolator  Form:  u^P^  ■  up _y  ♦  6  4  un  y  *  yfi2un1 

which  Is  independent  of  the  computational  path.  It  can  be  proved  that,  when  the 
trapezoidal  rule  is  used,  the  differential  partitioning  procedure  (9)  is  stable 
provided  the  extrapolator  is  chosen  in  the  form 


h2  - 

♦  h  Vi  ♦  r  Vi 


(in 


in  which  u  and  u  are  computed  by  di fferent iation  formulas  regardless  of  the 
the  implementation  paths  adopted.  The  above  extrapolator  for  the  differential 
partitioning  procedure  (9)  is  identical  to  that  of  path  ( C2 )  for  the  algebraic 
partitioning  procedure  (4)  (see  Table  4).  Hence,  it  appears  that  computational 
advantages  of  the  (JO),  (CO’),  (CO)  and  (Cl)  computational  paths  are  negated  in 
the  differential  partitioning,  whereas  they  are  preserved  in  the  algebraic  par¬ 
titioning. 
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PRESERVATION  OF  RIGID  BODY  MOTIONS 

The  spatial  convergence  of  nonconforming  finite  elements  has  been  tra¬ 
ditionally  tested  with  Irons'  patch  test  [9],  which  among  other  things  veri¬ 
fies  whether  rigid  body  motions  are  accurately  represented.  In  a  similar 
vein,  for  partitioned  solution  procedures  It  Is  Important  to  check  which  com¬ 
bination  of  extrapolator  and  computational  paths  preserves  rigid-body  motions. 
Note  that  If  the  temporal  order  of  accuracy  of  the  Implicit  Integration 
formula  is  two.  so  Is  that  of  partitioned  difference  equation  (4a).  Thus, 
one  is  tempted  to  believe  that  procedures  of  formally  second-order  accuracy 
should  trace  constant  velocity  and  constant-acceleration  motions  as  they 
vary  linearly  and  quadratical ly  In  time,  respectively.  The  formal 
order  of  accuracy  holds  true,  however,  only  for  solution  compo¬ 
nents  that  are  not  on  the  partitioned  boundaries  (see  Table  1).  Hence,  the 
partition  (4a)  does  not  In  general  preserve  the  two  practically  Important 
rigid-body  motions.  A  notable  exception  to  this  statement  occurs  for  those 
partitions  whose  matrix  admits  rigid  body  motions.  These  effects  will 
now  be  examined  assuming  that  the  trapezoidal  Integration  rule  Is  used. 

Rigid  body  motions  u  of  a  flexible  structure  satisfy 

Ku  •  0  (12) 

Introducing  the  partition  (4e): 

ll  «  *  12“  03) 

Now  equation  (13)  Implies  that 

(ft  ♦  «2  &i)  «  *  (H  -  «2  *2)  “  04) 
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If  computational  path  (O')  Is  followed,  the  following  difference 
equation  holds  (cf.  Table  2): 

<s“2 «,)  a„  •  «2  in  •  s<h„.,  ♦  a,.) 


2  P 
•  6  12  -n-1 


For  a  constant-acceleration  rigid-body  motion,  we  must  have 


05) 


H  ii  k  •  f  ,  k  •  0,  1  . . 


06) 


Inserting  (16)  Into  (15)  and  accounting  for  (14),  It  follows  that 
exactness  of  representation  demands  that  either  jC^  un  =  0  or  that 
the  following  predictor  be  used 


P 

*n 


H  n-1  *  26  -  n-1  4  26*  iin-1 


(17) 


This  Is  the  second  predictor  formula  given  for  path  (O')  In  Table  3. 

Remark  J.  If  one  Insists  only  on  exact  preservation  of  constant  velocity 
rigid  body  motions,  for  which  ufl  0  ,  the  following  predictors  can  be  used: 


P 

u 

-n 


“n-1  *  24  -  n-1 


(18) 


-n 


u.*2«0.+6  u, 
-n-1  -n-1  —  n-i 


Remark^J.  Both  the  element-by-element  and  node-by  node  partitions  preserve 
any  rigid-body  motions  because  these  two  partitions  satisfy  the  rigid-body 
conserving  Identity 


p 

For  these  partitions,  the  last-solution  extrapolator  u  ■  u  .  is  optl- 

—  n  —  n- 1 

mally  stable  for  path  (O')  and  preserves  rigid-body  motions.  For  other  com¬ 
putational  paths,  the  equivalent  predictors  listed  in  Table  ?  retain  these 
properties. 

Remark  3.  The  last-solution  extrapolator,  although  stable,  distorts  constant 
velocity  and  constant-acceleration  rigid  body  motions  for  the  DOF-by-DOF  and 
staggered  partitions.  This  effect  is  illustrated  in  Figure  1  for  a  two 
degree- of- f reedom  problem.  This  distortion  can  only  be  eliminated  by  iterating 
at  each  time  step  to  convergence. 


Remark  4.  The  accurate  extrapolator  (17)  gives  rise  to  numerical  instability 
for  paths  other  than  (O').  Thus,  stability  and  accuracy  (as  far  as  preserva¬ 
tion  of  rigid  body  motions  is  concerned)  are  not  generally  equivalent  attri¬ 
butes  for  partitioned  Integration. 


The  main  findings  of  this  study  are  sunmarired  in  Table  5.  Clearly 
path  (O')  is  the  choice  when  accurate  representation  of  constant-accelerat ion 
rigid-body  motions  is  Important,  and  a  DOF-by-DOF  partition  is  used.  The 
price  paid  for  this  formulation  is  the  need  for  additional  calculations,  such 
as  solving  for  accelerations,  in  the  advancing  step  [8],  If  the  accuracy  re¬ 
quirement  is  relaxed  to  preservation  of  constant-velocity  rigid-body  motions, 
then  either  path  (1)  or  (?)  is  acceptable  for  those  partitions  (because  pre¬ 
dictors  (18a)  and  (18b)  are  stable  for  those  paths,  cf.  Table  ?),  and  the 
computational  effort  Is  accordingly  reduced.  Finally,  path  (0)  distorts 
rigid  body  motions  other  than  constant  displacement  for  the  DOF-by-DOF  and 
staggered  partitions. 

So  far  we  have  focused  our  attention  on  the  construction  of  stable  and 
accurate  extrapolators  and  on  conditions  for  preserving  rigid-body  motions. 

The  remainder  of  the  paper  will  be  devoted  to  the  examination  of  the  “Incre¬ 
mental"  effect  of  partition  on  algorithmic  accuracy. 


DISPLACEMENT 


EFFECT  OF  PARTITIONING  ON  ALGORITHMIC  ACCURACY 


The  conventional  Fourier  technique  for  assessing  accuracy  of  Integra¬ 
tion  formulas  for  second-order  systems  proceeds  as  follows.  First,  the  homo¬ 
geneous  difference  equations  are  decoupled  by  projecting  them  on  normal  co¬ 
ordinates.  Second,  the  frequency  distortion  and  numerical  damping  of  an  un¬ 
coupled  difference  equation  subjected  to  harmonic  Input  of  frequency  w  are 
presented  as  functions  of  the  normalized  sampling  frequency  uh.  This  pro¬ 
cedure  Is  not  Inmedlately  applicable  to  the  partitioned  difference  equations 
(4b),  however,  because  these  are  not  generally  diagonalized  by  the  normal 
modes  of  the  semidiscrete  equations  (1). 

In  this  paper  a  limit  differential  equation  approach  will  be  utilized 
to  assess  the  frequency- dependent  accuracy  of  various  partitioned  procedures, 
assuming  that  the  solutions  can  be  expanded  In  Taylor  series  up  to  the 
order  of  the  Integration  formula.  This  approach  has  been  extensively  used 
for  the  evaluation  of  artificial  viscosity  [10,11].  It  will  be  shown  that 
the  limit  differential  equations  for  the  partitioned  difference  equations 
display  a  frequency  distortion  effect  as  primary  error  source  due  to  par¬ 
titioning. 

As  an  Illustration  when  the  trapezoidal  formula  Is  used  In  conjunc¬ 
tion  with  path  (CO*)  and  the  previous  solution  vector  u  n Is  used  as 
the  extrapolator,  the  following  limit  differential  equations  can  be  derived 
(see,  for  example,  Roache  [11]  for  details), 

(H  -  7<2  il  4  I®2  H  4  £,  M  "  0(*3)  0*) 

Note  that  the  mass  matrix  Is  considerably  modified  whereas  the  stiffness 
matrix  Is  not,  and  also  the  absence  of  terms  associated  with  the  velocity 
vector  u  .  We  therefore  conclude  that  the  primary  additional  algorithmic 
error  caused  by  partitioning  Is  manifested  In  the  frequency  distortion, 


which  has  to  be  evaluated  from  its  frequency  equation 


det  | (M 


2  2 

56 


K  i 
**  I 


«  K?)  n 


2 

6  K  |  ■  0  , 


a  ■  u  6 


(20) 


The  two  limit  cases  for(20)  are:  *  0  (fully  implicit  integration)  and 

X->  ■  K  (fully  explicit  Integration).  These  correspond  In  turn  to  the 
standard  trapezoidal  rule  and  central  difference  formula,  respectively.  It 
is  well  known  that  these  two  formulas  Introduce  no  numerical  damping.  More¬ 
over,  the  central  difference  formula  shortens  the  period  whereas  the  trape¬ 
zoidal  rule  elongates  It  by  an  amount  roughly  twice  as  big;  this  can  be 
imnedlately  deduced  from  (20).  That  the  frequency  distortion  Is  also  the 
primary  algorithmic  error  due  to  partitioning  can  also  be  shown  to  be  true 
for  other,  numerically  damped.  Integration  formulas.  (The  demonstration  re¬ 
lies  on  the  fact  that  nonzero  odd-derivative  residual  terms  In  the  Taylor 
series  occur  beyond  the  order  of  the  Integration  formula  truncation  error). 

As  noted  previously,  the  matrix  equation  (19)  Is  not  generally  dlagona- 
lizable  by  the  natural  modes  of  the  original  equations  of  motion.  To  assess 
the  accuracy  performance  of  various  partitions  we  can  exploit,  however,  the 
fact  that  the  frequency  error  due  to  partitioning  emanates  from  partition 
boundaries.  This  enables  us  to  construct  a  two  degree- of- freedom  model  system 
and  to  assess  the  fequency  distortion  due  to  partitioning  from  that  of  the 
model  system.  The  governing  equations  of  this  system  are 


1  ♦  a 

0 

*• 

u  ♦ 

1  ♦  - 
a 

-  1 

0 

1 

-  1 

1 

with  appropriate  Initial  conditions.  Equations  (21)  can  be  viewed  as  the  axial 
equations  of  motion  for  a  fixed-free  bar  consisting  of  two  elements  whose 
length  ratio  Is  1  :<*  . 

For  the  above  two-d.o.f.  system  the  partition-caused  artificial  fre¬ 
quency  distortion  can  be  evaluated  by  equation  (20).  Alternatively,  the 
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distortion  can  be  computed  exactly  from  the  characteristic  equation  for  the 
partitioned  difference  equations  (see  Table  3).  As  the  calculations  by  the 
two  approaches  have  shown  no  discernible  difference,  we  present  the  fre¬ 
quency  distortion  by  the  latter  approach  in  the  sequel. 

Fiqures  ?a  and  ?b  show  the  frequency  distortion  of  the  low  and  high 
frequency  components  for  the  equal  length  case  n  ■  1 .  for  this 
case  thp  node-by-node  imp) ic it-expl ic  i t  procedure  is  the  most  accurate,  fol¬ 
lowed  by  the  element-by-element ,  the  DOF-by-DOF.  and  the  stagqered  procedures. 
As  the  element  lenqth  ratio  (roughly  square  of  the  frequency  ratio)  is  de¬ 
creased.  however,  an  intermediate  value  of  o  is  reached  at  which  the 
element-by-element  procedure  becomes  more  accurate  than  the  node-by-node 
procedure  (see  Fig.  3).  If  the  length  ratio  is  further  decreased,  the  low 
(high)  frequency  error  approaches  that  of  fully  explicit  (implicit)  formula, 
respectively,  as  can  be  shown  from  Figure  4.  The  increase  in  the 
frequency  ratio  weakens  the  coupling  effect  between  the  two 

From  the  foregoing  accuracy  analysis  of  the  model  system  it  can  be 
seen  that  accuracy  of  solution  components  far  from  the  partition  boundaries 
is  larqely  controlled  by  that  of  the  integration  formulas  used  For  solution 
components  at  or  near  the  partition  boundaries,  the  element-by-element  and 
the  node-by-node  partitions  enjoy  a  somewhat  higher  accuracy  than  either  the 
explicit  or  the  implicit  formula  itself  in  the  sense  that  the  frequency  dis¬ 
tortions  are  less  than  those  of  the  fully  explicit  (implicit)  integration 
formulas  (see  Fiqs.  2  and  3).  On  the  other  hand,  for  the  DOF-by-DOF  and  the 
Staggered  partitions  the  accuracy  of  the  boundary  solution  components  is 
inferior  to  that  of  the  fully  explicit  (implicit)  formulas. 
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CONCLUSIONS 


The  material  covered  herein  complements  an  earlier  paper  [1],  which  was 
devoted  to  the  qeneral  categorization  of  partitioned  integration  procedures. 
The  present  paper  expands  on  stability  and  accuracy  as  properties  strongly 
related  to  the  choice  of  extrapolation  formula  and  computation  path;  these 
two  aspects  were  only  briefly  commented  upon  in  [1].  The  information  pre¬ 
sented  here  may  be  used  as  an  initial  quide  in  method  selection  process.  It 
whould  be  noted,  however,  that  the  selection  of  one  partitioned  procedure 
over  another  may  be  based  on  considerations  other  than  those  discussed  here, 
for  example,  the  availability  of  computational  tools  for  treating  certain 
subsystems  of  the  overall  problem  (e.g.,  finite  element  structural  analyzers, 
boundary-integral  analysers  for  Infinite  domains,  etc.)  may  naturally  dictate 
the  use  of  certain  partitions. 

We  have  sunmarl zed  In  Table  6  qeneral  algorithmic  characteristics  of 
five  partitioned  solution  procedures,  from  an  unconditional  stability  view¬ 
point  the  Implicit-implicit  procedures  (staqqered  and  element-by-element 
partitions)  would  be  preferred.  This  would  be  the  case  if  the  envisioned 
time  stepsize  is  larqe  compared  to  the  shortest  characteristic  times  of 
each  subsystem. 

Implici t-implicit  partitions  are  desirable  if  subsystems  display  widely 
different  response  characteristics .  They  offer  the  advantage  of  computa¬ 
tional  simplicity  for  the  expl icitly-treated  subsystems.  From  an  applica¬ 
tions  viewpoint  the  element- by-element  partition  is  natural  for  finite  element 
discretizations,  but  can  become  ci/nbersome  for  extensively  connected  parti¬ 
tions  (e.g.,  surface  of  volume -coup led  finite  element  models).  The  node-by¬ 
node  partition  is  estimated  to  require  additional  implementation  effort,  as 
"boundary  nodes'  must  be  appropriately  labeled.  It  becomes  attractive,  how¬ 
ever,  when  the  number  of  boundary  unknowns  is  considerable.  The  DOF-by-DOF 


partition,  althouqh  less  accurate  than  the  previous  ones,  may  be  advantageous 
when  the  explicit-implicit  partition  is  to  occur  within  each  element  or  node 
point  of  a  finite  element  or  finite  difference  model.  An  example  would  be 
the  computational  separation  of  translational  and  rotational  degrees  of  free¬ 
dom  in  lumped-mass  finite  element  models. 

We  now  simmarize  our  main  findings. 

For  a  given  integration  formula,  stability  and  accuracy  of  partitioned 
integration  procedures  are  affected  by  the  interaction  of  the  computat ional 
path  and  selected  extrapolator.  The  user  is  free  to  select  partition-type . 
computational  path,  and  extrapolator  from  many  possible  combinations.  Such 
versatility  is  an  important  attribute  of  the  present  formualtion:  implicit 
time-discretization,  algebraic  partitioning,  and  extrapolation.  It  is  not 
shared  by  the  "differential  partitioning"  formulation:  partitioning  at 
differential  equation  level,  extrapolation,  and  time  discretization,  which 
was  used  in  [2-6]. 

Possible  distortion  of  rigid-body  motions  due  to  partitioning  can  be 
prevented  for  most  partitions  if  the  computation  path  (CO*)  in  conjunction 
with  a  proper  extrapolator  are  adopted.  This  is  not  the  case  when  the  dif¬ 
ferential  partitioning  approach  is  used  except  in  the  case  of  element-by¬ 
element  and  node-by-node  partitions. 

The  primary  error  due  to  partitioning  is  manifested  in  frequency  dis¬ 
tortion  for  solution  components  adjacent  to  a  partition  boundary.  Partition- 
caused  nianerical  damping  is  of  secondary  importance.  In  particular,  for 
integration  formulas  which  possess  no  numerical  damping,  such  as  the  trape¬ 
zoidal  and  the  central  difference  formulas,  partitioning  Introduces  no  addi¬ 
tional  numerical  damping  if  the  extrapolator  is  judiciously  selected. 
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